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PURPOSE 


In dealing with large two-dimensional arrays (such as those representing 
digitized pictures), one difficult problem is to determine quickly and cheaply 
the result of an operation. For example , during the development of a program 
to filter digital data, the logic of how to handle the edges, loss of sync, 
dropped lines and so on is complicated enough to cause major debugging prob- 
lems. A simple , even well known,: program such as one to rectify or register 
data must be tested and debugged on the particular format of the data avail- 
able. The whole process of developing software is slowed when a special 
display device must be used to view an intermediate stage of a program. It 
has long been recognized that a natural way of obtaining a quick (if crude) 
view of a large array is to use the high-speed printer of the computer with 
different symbols representing ranges in the data. 

The problem attacked by program OCM (One Channel Map) has a more critical 
feature. In order to obtain training data for a classification algorithm, 
the exact (digital) coordinates of the boundaries of individual fields must 
be known. Ibr this reason, even high quality hard photographic information 
may be useless. On the other hand, average quality programs (to produce 
maps on the printer) fail to display the data in sufficient detail to allow 
the user to locate (on the printout) fields corresponding to ground-truth. 
Program OCM fills the need for better quality reproduction of digital pictures 
on the high-speed printer. 

METHOD 

Program OCM achieves the objective by printing an appropriate selection 
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of characters eight lines per inch on blank (unlined) paper. The translation 
of the data is performed by a subroutine written by Mr. Dale Ruspino (at 
Texas A&M) which uses effectively the IBM System/360 "translate under mask" 
instruction to translate long strings in place, essentially instantaneously. 
Experimentally it was learned that, what ever the choice of symbols, only five 
levels of gray are unambiguously detected with a single print, and as the 
ribbon wears this number becomes four. With a new ribbon and overprinting 
once, with an additional translate to get the overprinted character, eight 
levels are detectable. 

In order to use the translate subroutine, translate tables must be set 
up. The minimum and maximum significant values in the data must be specified. 
The program divides the data spread into eight equal (as possible) partitions, 
sets all lower values in the first partition, all higher in the eighth, and 

assigns symbols (suggested symbols are given in USAGE below) to each range. 

‘ ^ 

(The value 0 is assigned a special symbol.) The data is assumed to be a 
positive, in that low values represent dark areas; thus the first three 
symbols are overprinted to make them darker. 

Program OCM then initializes a variable which breaks the data into one 
page size parts of not more than 126 columns each. (If more than this many 
columns are required, the program will return to this point, continuing until 
all data is printed.) 

The program then initialized the associated variable IREF for direct 
access device 1 (set up by a DEFINE FILE statement and appropriate JCL) , 
enters a loop to read, find the next row, translate and print and translate 
and overprint the row just read, until all rows are finished. If any columns 
remain, the program returns for more printing. 
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The translate program Is so fast that no difference can be detected 


whether It Is present or absent by looking at the run time. 
DESCRIPTION OF PROGRAM VARIABLES 


DUM(4) 

LOGICAL*! 

Dummy vector equivalenced to iT - used In definition 
of translate table TS. 

I» J»L 


Generic DO lndicees 

ICilCP 


First and last column index currently being printed 

ICE, ICS 


Last and first column index (read as data) to be 
printed 

ICHAN 


Channel number (read as data) 

IRE, IRS 


Last and first row index (read as data) 

IREF 


Associated variable for unit 1 

IT 


Integer used as an index by filling the last byte 
via DUM and an equivalence statement. 

IW 


Write parameter used in setting up translate table TI 

MAX, MIN 


Maximum and minimum expected significant data values 
(read as data) 

R(500) 

LOBICAL*l 

Buffer into which each row is read. 

SYMBOL (12) 

LOGICAL* 1 

Symbol table: first for 0, next eight first break- 
down (i.e. for TR)next three overprinted over 2,3,4 
(via TS) 

TR(256) 

LOGICAL*! 

Translate tables for subroutine 

TS (256) 

LOGICAL*! 

HP 05 

X8 


(Max-Min) / 8 . j used to accurately spread the data 
into eight : parts . 



USAGE 


Detailed deck set up and JCL: TAMU 360/65 with HASP 

//NAME J0B (standard jobcard) 

HASP CARDS: 

/*CLASS A (0-110K) 

/*F0RMS 14101010 (blank paper , 1 part, 8 lines/inch setup) 

/*LINES/PAGE 0 (user controls spacing to new page) 

/*R0UTE PR1NTER1 (selects a generally better operated printer) 

EXEC CARD; 

// EXEC F0RTG,REGI0N= 110K (invokes FORTG cataloged procedure) 

DD CARDS: 

//SYSLIB DD 

// DD 

// DD DSNAME=USER.T4O5.EE.PR0BUB(HPO5),DISP=SHR 

(this JCL allows FORTRANG to retain its SYSUB and adds HEX 
program HP05, resident on the described private disk) 

//FT01F001 DD UNIT- 2314,V0L=SER=JACKO1,DSN=HILLC0, 

// SPACE- (500, 12000) ,DISP=(0LD,KEEP) ,DCB= (RECFM-F ,DS0RG=DA) 

(this DD card makes F0RTRAN UNIT 1 a direct access device) 

/ /S0URCE DD * (source program next) 

PROGRAM 0CM DECK 

//SYSIN DD * (data next) 

DATA (described next) 

/* (end of file) 
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DATA DESCRIPTION: 


SYMBOL FORMAT 12A1 

I CHAN , IRS , IRE, ICS , ICE FORMAT 513 

MIN .MAX FORMAT 213 

repeat for as many pictures as may be desired; each picture thus takes two 
data cards. 

Subroutines required: 

HP05 Assembly language program to perform translation. 

C0LS Program to print column headings, which requires 

T0DEC Program which furnishes the digital (base 10) representation 

of the column Index. 

Restrictions : 

This program is highly machine dependent, even some what installation 
dependent. Clearly, however, the idea can be implemented, for instance on 
the 1108, using JjTRAN and a translate instruction subroutine PRETTY (documented 
in the Lunar Orbiter Project) written in SLEWTH. The logic is trivial in any 
case . 

The use of this program requires the data be written on a disk one line 
of 500 bytes each (125 words) per record unformatted. The data is expected 
to be lined up somewhat differently than it was on the tapes; there are 
12000 records with the J-th row of channel I having index I + 24*(J-1). 

That is, row 1 channel 1 through row 1 channel 24 then row 2 channel 1 and 
so on. The value 0 represents no data and is printed differently. 
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Suggestions: 

SYMBOL: A convenient selection of symbols to print is 
U (for zero) 

M overprinted with $ 

0 overprinted with + 

0 overprinted with - 

* 

+ 

blank 

It has been determined that appropriate values for MAX and MIN for the display 
of field boundaries are about the top and bottom 7 to 10 percent on a cumulative 
distribution of all data values. Quite possibly other purposes will require 
different selections of MAX and MIN. 

Run time: 

Approximately 0.05 seconds per 126 character row desired. Thus, a run 
of 12 channels, cols 50 to 427 (requiring 3 126 col pages to print), rows 80 
to 458 should take about 682.20 seconds; the actual run time was 689.40 seconds, 
including overhead charges. This is a lot of printing for which the overhead 
is high. 

Lines output: 

Specify two lines output for each line actually printed, plus an allow- 
ance for column headings and program and JCL listing. 



REFERENCE INFORMATION INCLUDED: 


General flow chart 
Detailed flow chart 
Documentation of COLS. 
Comment on HP05 


Sample run 



/ REBP SYMBOLS 
REBP CHANNEL NUMBER/ 
START 4 END ROW 
5TRRT 4 END COL . 
RCRP MIN BNP MBX / 



' PRINT CHANNEL NUMBER/ 
CDL NUMBERS ON THIS PAGE 
PRINT COL HEADINGS / 



/TRANSLATE IN PLACE 
PR I NT 

TRANSLATE IN PLACE, 
OVERPRINT / 


10RE ROWS'; 


GO TO MORE ROW 


C> A. | u 




O C/'t ■ 
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SUBROUTINE DOCUMENTATION : SUBROUTINE COLS 

PURPOSE: Subroutine COLS prints column headings at the top of a page written 

vertically with leading zeros suppressed. It is used to greatly simplify 
the finding of the coordinates of field boundaries. 

METHOD: In order to both speed the operation and simplify the code two arrays 

are employed. The full word Integer vector FAST is equivalenced to a 
L0GICAL*1 array PAGE. Initially FAST is filled with all blanks (40 Hex). 

Then a loop is entered which fills FAST with the digits of the column headings 
to be pripted in the appropriate positions. Logic is included which suppresses 
leading zeros. The conversion to base ten is performed by subroutine T0DEC 
(which is self-explanatory). Then FAST is printed in one 1/0 statement. 

USAGE: 

Calling sequence: Call C0LS(I,J,N) 

where I is the index of the first column to be printed 

4 

J is the index of the last 

N is the printer position corresponding to the first column 

Restrictions: Obviously, 0SI£JSI+ 132-N. A diagnostic to check this 

is included in the program. Also, T0DEC requires J £ 999. 

REFERENCE INFORMATION: Detailed flow chart, C0LS and T0DEC. 


IS 




cu s 









WRITCCG/ I I ) 
ERST 


(return:) 




WR J TECG/2> / 
1/ J, N / 

T , 

(return) 


EORMRT CODES** 

I l«C ' I '/33AH/' '/ 33RH/ ' «,33AH/> 

2»( ' UNRCCEPTR0LE INPUT TO C0U5! '/ 33 10/' I » ) 

PROGRAM VRRIRBLE5 D AND BL ARE INITIAL 1 2ED IN TYPE 
DECLARATION STATEMENTS. VARIABLES PAGE AND EAST ARE 
EflUlVALENCED. 


1 7 


P. i *. t 


-f ri r( f 


C / 




FORMAT STATEMENT REFERENCED : 

2-C UNACCEPTABLE INPUT TO TODEC — FATAL ERROR’/ IJ0> 


iS 

// r .-.y r/rc.Y 7 C<r>r C 



SUBROUTINE COMMENT : SUBROUTINE HP 05 

Authors: Dale Ruspino and Bruce Marion, July 1974. 

Calling Sequence: CALL HP05(I,R,TR) 

where 

I is the length in bytes of that portion of the vector R to be 
translated; 

R is a vector of arbitrary length (within the maximum dimension 
restriction of FORTRAN) to be translated in place; 

TR is the vector of length 256 (bytes) containing the translate table. 

Method: Using the 360 assembly language instruction TM, the vector R is broken 

down into 255 length pieces, translated (ie the bytes are replaced by others 
according to what is contained in TR) until all I bytes have been processed. 

For example, suppose 

I is 2, R is (in Hex) 

01 0D and TR starts out 

CC 17 DD 40 40 40 40 40 40 40 40 CC 22 13 ... 

After the call to HP05, R is (in Hex) 1713. HP05 attains its speed by the 
use of this parallel translate instruction and by not saving registers which 
are not disturbed and not checking the input for errors. Indeed, HP05 has 
proved to be at least as fast as a dummy program^(in FORTRAN) 

SUBROUTINE HP05(I,R,TR) 

LOGICAL*! R(l) ,TR(1) 

RETURN 

END 

used while debugging the main program. 

Reference Information: Listing 
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PURPOSE 


In drawing one channel maps on the high-speed printer, one needs to set 
the maximum and minimum ejected significant data values. It turns out that 
for the 24 channels we have to work with these values vary wildly; for example 
channel 22 has its 5 to 95 percent cumulative distribution values 27 and 
62, whereas channel 24 (on the same pass) has the same values 20 and 31. 
Program HIST06 accumulates, prints and punches complete histograms of all 
24 channels. The card output is available in easy to read form for further 
analysis. 

METHOD 

The program collects data stored on a direct access device and prepares 
histograms. Since the information is packed by rows, it must be unpacked be- 
fore it can be used as integer information. This function is performed in 
FORTRAN by filling the fourth byte of a dummy logical vector (L0GICAL*1) 
which is alllgned with a full word integer using an EQUIVALENCE statement. 


DESCRIPTION OF PROGRAM VARIABLES : 

H(128) HJTEGER*4 Used to accumulate histograms 


I,J , JJ ,JP 
IF 


Generic DO loop parameters 

Used in FIND instruction on direct access 
device unit 1 


IR 

L(4> L0GICAL*1 
M INTEGER*4 
LOGICAL*! 


Associated variable for unit 1 
L and M are assigned the same storage 
location; used to unpack information 
Read buffer-contains row most recently read. 


R(500) 



USAGE 


Detailed deck setup and JCL:TAMU 360/65 with HASP 
//NAME JOB (standard jobcard) 

/*CLASS A 

// EXEC F0RTG,REGI0N=11O K 
//FT01F001 DD UNIT-2314, V0L-SER-JACKO1, 

// DSN-HILLC0 , SPACE- (500, 12000) ,DISP- (0LD , KEEP ) , 

// DCB- (RE CFM-F , D S 0RG-DA) 

//S0URCE DD * 

HIST0G PR0GRAM DECK 
//SYSIN DD * 

/* j 

(For more detailed comments on JCL, see the documentation of program OCM.) 
Description of output: The main output is the card output. 

(The same information is printed, however.) Each channel produces 15 cards, 
as follows: 

CHANNEL-SEQUENCE NUMBER (15) 

ONE SPACE (IX) 

H(1),...,H(9) (918) 

repeat CHANNEL-SEQUENCE NUMBER, IX 

H(10) , . . ,H(18) (918) 

The fifteenth (last) card for each channel contains the CHANNEL-SEQUENCE 

NUMBER, IX and H(127) , H(128) . 

Restrictions: The data value 0 is Ignored altogether. Values over 127 



(none were found) are accumulated In H(128) . 

Run time: Slightly less than six minutes are required for all 24 channels; 

this amounts to about 60 jis for each value unpacked and accumulated, counting 
all overhead. The overhead associated with 12000 reads Is, however, consider- 
able; probably less than a third of the time shown on the printout la spent 
doing calculations. 

REFERENCE INFORMATION INCLUDED 


Detailed flow chart 


Sample run 
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PURPOSE 


Program HISTOG accumulates raw data (histograms) from the 24 channels 
stored on the disk; what is needed, however, are estimates of the minimum and 
maximum significant data values present as input to program OCM. Program 
PRINTUM calculates the cumulative distributions and also prints histograms 
(on the printer). Only the 5 to 95 percent distributions are printed. 

METHOD 

Program PRINTUM reads the channel number and histogram vector H, blanks 
the LOGICAL *1 array PAGE (126,127), determines the maximum value M present 
in the data, fills page with out to the properly scaled index H(I)*126/M, 

and then prints the histogram. The total number of values accumulated is 
counted and a cumulative distribution is taken of H (starting with value 1) ; 
those percentages lying between 5 and 95 percent are printed. The program 
returns for more data until an E0F Is encountered. 

DESCRIPTION OF PROGRAM VARIABLES 

Note: all are INTEGER*4 unless otherwise indicated. 

A DO loop index 

B L0GICAL*1 Blank 

C Channel number; read along with H 

H(127) Vector - the raw histogram 

I Generic DO loop index; then used to print descending 

J Generic DO loop index and parameter 

JJ, JP Generic Do loop index and parameter 

M Maximum value found in H 
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PAGE (126, 127) 

L0G1CAL*1 Array allowing print of histograms 

R 

Percent printed 

S 

DO loop index; then the grand total of values in H. 

ST 

LOGICAL* 1 '*’ 

T 

Cumulative sums of H 

V 

DO loop index. 

USAGE 


Detailed deck setup and JCL: TAMO 360/65 with HASP 

//J0BNAME 

J0B (standard job card) 

/*CLASS 

A 

/★LINES /PAGE 

0 

// EXEC 

F0RTG , REGI0N** 110K 

//S0URCE 

DD * 

PR1NTUM PR0GRAM DECK 

//SYSIN 

DD * 


DATA- THE CARD 0UTPUT 0F HIST0G 


/* 

Data description: 

The card output of HIST0GJ specifically, each card contains the channel 
number (13) a sequence number which is skipped (3X) and data for reading 
into H (918). 

Lines Output: 

Approximately 4000 lines for 24 channels. 



Execution time 


Approximately 3.7 seconds per channel plus overhead; a total of 1.63 
minutes for 24 channels. 

REFERENCE INFORMATION INCLUDED : 

Detailed flow chart 
Sample run 
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